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We present in detail the implementation of the Blaizot-Mendez-Wschebor (BMW) approximation 
scheme of the nonperturbative renormalization group, which allows for the computation of the full 
momentum dependence of correlation functions. We discuss its signification and its relation with 
other schemes, in particular the derivative expansion. Quantitative results are presented for the 
testground of scalar 0{N) theories. Besides critical exponents which are zero-momentum quantities, 
we compute in three dimensions in the whole momentum range the two-point function at criticality 
and, in the high temperature phase, the universal structure factor. In all cases, we find very good 
agreement with the best existing results. 



I. INTRODUCTION 

The exact or non perturbative renormaliza- 
tion group (NPRG), as formulated in the sem- 
inal work of Wetterich [ij , leads to an exact flow 
equation for an effective action (see also Refs. 
[23 for the original formulation of the NPRG). 
This equation cannot be solved in general, but of- 
fers the possibility of developing approximation 
schemes qualitatively different from those based 
on perturbation theory, allowing us in particular 
to tackle non-perturbative problems. 

The "derivative expansion" (DE) is such a 
scheme: based on an expansion of the running ef- 
fective action in terms of gradients of the fields, 
it has been applied successfully to a variety of 
physical problems, in condensed matter, parti- 
cle physics or statistical mechanics (see e.g. [4]). 
The DE scheme allows us to calculate not only 
universal, but also non-universal quantities de- 
fined at vanishing momenta, such as critical ex- 
ponents and phase diagrams (see e.g. 04^)- 
However, it docs not give access to the full 
momentum dependence of correlation functions, 
something desirable in many situations. 

The present paper deals with another scheme, 
namely with the strategy proposed by Blaizot, 
Mendez-Galain and Wschebor (BMW) in [T^, Igj 



which is reminiscent of earlier attempts by Parola 
and Reatto The BMW scheme makes use of 
the cutoff function Rk{q), typical of NPRG stud- 
ies, which renders all vertex functions smooth, 
and insures that the fiow at scale k involves the 
integration of fluctuations with momenta q at 
most of order k. At order s, the approximation 
consists in setting the internal momentum q to 
zero in the vertex functions of order larger than 
s, leaving a closed set of flow equations for the 
first s vertex functions. It was shown in [7] that 
the BMW method encompasses any perturbative 
results, provided it is pushed to high-enough or- 
ders s. For the particular s — 2 case studied 
in the following it is one-loop exact for the two- 
point function. It has also been shown in Q that 
in the N ^ 00 limit of the 0{N) scalar models 
the BMW scheme becomes exact and allows for 
the computation of any correlation functions. 

The BMW approximation scheme has been ap- 
plied to 0{N) models with success, in simplified 
versions involving either expansions in the fields 
[lo| or an approximated propagator fll|. In this 
paper, we provide a full account of the practi- 
cal implementation of the method without fur- 
ther simplifications. We also detail and extend 
the first results obtained recently at order s — 2 
0]. In particular, we extract the universal scal- 
ing function governing the critical region. At the 
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quantitative level, we find very good agreement 
with the best existing results. 

The remainder of this paper is organized as fol- 
lows: after a general presentation of the NPRG 
framework (Sect|ll|, we detail the BMW approx- 
imation scheme in Section IIIH and show how to 
implement it in practice in Section IIVI In Sec- 
tion |Vl we report on the results obtained at crit- 
icality (critical exponent values, shape of two- 
point function, etc.), while in Section IVll we com- 
pare the universal scaling function obtained in 
the whole critical region to existing results. In 
Section we show how the BMW results shed 
new light on the derivative expansion approxi- 
mation scheme, and we discuss in particular its 
validity domain. Conclusions can be found in 
Section I Villi while the appendices gather tech- 
nical material. 



II. THE NPRG FRAMEWORK 

For the sake of simplicity, in the following, we 
only discuss a scalar field theory in the Ising uni- 
versality class, and defer to Appendix ID] the pre- 
sentation of the equations that hold for general 
0(A^)-symmetric models. 

We consider the usual partition function 



Dip e 



with the classical action 
1 
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■ if 



4! 
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In Eq. ([T]), j{x) is an external source and /^.^V 
is a shorthand for J d'^x j{x)(p{x). 

The NPRG strategy is to build a family of the- 
ories indexed by a momentum scale parameter fc, 
such that fiuctuations are smoothly taken into ac- 
count as k is lowered from the microscopic scale 
A down to In practice, this is 

achieved by adding to the original Euclidean ac- 
tion S a fc-dependent quadratic (mass-like) term 
of the form 



with 



d\ 



(3) 



so that the partition function at scale k reads 



Zk[j] = J V^e-^'-^s^+f^^^ . (4) 

The cut-off function Rk{q) is chosen so that: i) 
it is of order k^ for q k, which effectively sup- 
presses the modes Lp{q k); ii) it vanishes for 



q ^ k, leaving the modes (p{q) with q k un- 
affected. Thus, when fc = A, Rk{q) is of order 
A^ for all g < A, and fluctuations are essentially 
frozen. On the other hand, when fc = 0, Rk{q) 
vanishes identically so that Z^^q — Z, and the 
original theory is recovered. The specific form of 
the cut-off function Rk{q) will be specified later. 

Following Wetterich [Ij, an effective action at 
scale fc, Tfc [(/)], is defined through the (slightly 
modified) Legendre transform 



Tkick] + logZfc[j] - /" j> - A5fe[0], 

J X 



(5) 



where 4){x) = 5\nZk[j]/ 5j{x). This effective ac- 
tion obeys the exact fiow equation [l| (up to a 
volume factor): 



dtVk 



dtRk{q) 



g [g,-(j;0]-f i?fe(q) 



(6) 



where dt = kdk and T]^'[q,—q](j)\ is the Fourier 
transform of the second functional derivative of 



T^(2)r 



6^Tk 



54){xi)5(t){x2) 



(7) 



Thus, (rf '[g, -q- cj)\ + Rk{q)y^ is the fuU propa- 
gator in the presence of the field (j){x). The initial 
condition of the fiow equation ([SJ is specified at 
the microscopic scale fc = A where fiuctuations 
are frozen by ASk, so that rfe=A[0] « S[cf)]. The 
effective action r[^] of the original scalar field 
theory is obtained as the solution of Eq. ([6]) for 
fc — 0, at which point Rk{q) vanishes identically. 

When (f) is constant, the functional Tk[(t>\ re- 
duces, to within a volume factor f2, to the effec- 
tive potential Vk{4>)'- 

Ffc 17 Vfc(0), constant. (8) 

The flow equation for Vk reads 

dtVk{p) = \ j dtRk{q)Gk{q,4>), (9) 



where 



GZ\q,4>)^V^;^\q,4>) + Rk{q), (10) 



(see Appendix lAl for notation). 

Let us now consider the flow of n-point func- 
tions. By taking two functional derivatives of 
Eq. ([U, letting (f) be constant, and Fourier trans- 
forming, one obtains the equation for the 2-point 
function: 
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a^rf (p, </)) = y dtRk {q)Gl {q, (l>r^'> {p, q, ~p - q, 0)Gfe {q + p, ^)rf {~p, p + q, -q, <f>) 



(11) 



The flow equations © and PT|) are the first 
equations of an infinite tower of coupled equa- 
tions for the n-point functions: typicafly the 

(n) 

equation for T^, involves all the vertex functions 
up to r("+^)_ Approximations and truncations 
are thus needed to obtain any practical result. 
The presence of a sufficiently-smooth cutoff 

function Rk{q), (i) insures that the r[."^'s re- 
main regular functions of the momenta and (ii) 
limits, through the term dtRk{q), the internal 
momentum q in equations such as Eq. (jlip . to 
q k. These key remarks allow for approxi- 
mations without equivalent in more traditional 
frameworks and thus constitute one of the speci- 
ficities of the NPRG approach. 

The approximation scheme most widely used is 
the derivative expansion, which is entirely based 
on the above remarks about the analyticity of the 
vertex functions. It amounts to formulating an 
ansatz for Tk[4>] as an expansion in the deriva- 
tives of the field. For instance, at order V^: 

(12) 

The fiow equation ([6]) then reduces to a set of 
two coupled, partial differential equations for the 
functions and Zk{(t))- The derivative ex- 

pansion scheme has produced, along the years, a 
wealth of remarkable results (see e.g. d, d, [Tsj ) 
but it does not allow to access the full-momentum 
dependence of the vertex functions, something in- 
herently possible within the BMW scheme. We 
discuss this point in Section [VIIl where we show 
how the BMW approach sheds new light on the 
derivative expansion. 

III. THE BMW APPROXIMATION 
SCHEME 

The BMW scheme at order s aims at pre- 

(s) 

serving the full momentum dependence of T\. 
and approximating the momentum dependence 
of rl'^^'' and v'f^'^'^ in the fiow equation of v\^^ 

M 

For uniform fields, the following formula: 

ri^^'Hfe},o,0) = a^ri^)(te},0) (la) 

where the index i runs between 1 and s, allows to 



reduce the order of the vertex functions as soon 
as one momentum is vanishing. 

The BMW approximation relies on this for- 
mula, together with the analyticity of the vertex 
functions and the fact that the internal momen- 
tum q in the flow equations such as Eq. pT|) is 
effectively limited io q < k. The BMW scheme 
at order s thus consists in: 

(i) neglecting the dependence on the internal 
momentum q of F^^"'^^'' and F^""*^^': 

^k^^\pi,---,Ps -q,q,^) 

f(^+i)(pi,...,P„O,0), (14) 

and similarly for r\^^^\pi, . . . ,ps,-q,q,(j))] 

(ii) using Eq. which allows us to express 
the approximated expressions ((Ti)) as derivatives 

is) 

of F^ with respect to 0, thereby closing the hi- 
erarchy of RG equations at the level of the flow 
equation for F^*''. 

Note that the substitution in Eq. (fT4)) is not ap- 
plied to the g-dependence already present in the 
bare n-point functions [7, 8]. Thus, for instance, 
at the lowest level of the approximation (s = 0, 
the local potential approximation discussed be- 
low in Sect. lIIlX| . one leaves untouched the bare 

dependence oiT^^\q). This ensures in partic- 
ular that the propagator is one- loop exact. 

The accuracy of the scheme depends on the 
rank s at which one operates the approximation. 
Obviously the implementation becomes increas- 
ingly complicated as s grows. We will show later 
that good results can be obtained with low order 
truncations, i.e., at the levels s — Q and s = 2. 
The corresponding approximations are discussed 
in the next subsections. 



A. s = 0: The local potential approximation 

The local potential approximation (LPA) is of- 
ten seen as the leading order of the DE approx- 
imation scheme [1, [ij. In this subsection, we 
show that it can be seen also as the zeroth order 
of the BMW scheme. 

The BMW approximation for s = 0, consists 
in neglecting the (nontrivial) g-dependence of the 
2-point function in the flow equation ^ of the 
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"zero-point" function, that is of the effective po- 
tential Vfe. That is, one substitutes 

rl^) {q, c^)^q' + r[^^ (0, 0) - + dlv, .(15) 

Note that the equahty in the equation above is 
a particular case of the general relation ((T3)) . By 
substituting Eq. ([T5|) in Eq. ^ , one gets the equa- 
tion for the potential in the form 

This is the flow equation for the potential ob- 
tained within the DE truncated at the LPA level. 
There, Eq.® is derived by computing the prop- 
agator from the ansatz 

rLPA[^] = 1 d'^x |1 iv^f + yfc(0)| (17) 

and inserting it Eq.®. 

Since it allows for the calculation of the entire 
effective potential, the LPA provides information 

(n) 

on an the "sat once but only for vanishing ex- 
ternal momenta: these functions are indeed those 
that are obtained by taking the derivatives of the 
effective potential, i.e., 

r^^\Or-- A<l^) = d^Vk. (18) 

Non trivial momentum dependence will appear at 
the next level of approximation, to be described 
in the next subsection. 



B. First order with full momentum 
dependence: s = 2 

The order s = 2 is the first order of the ap- 
proximation where a non-trivial momentum de- 
pendence is kept. The loop momentum q in the 
3 and 4-point functions in the right hand side 
of Eq. PT|) is neglected, and Eq. is applied. 

(2) 

The flow equation for Fj. {p,(j)) becomes then a 
closed equation 

(19) 

where we have introduced the notation: 

/„(0) = J„(p = O,0), 
J„(p,0)= / dtRkiq) Gk{p + qA)Gl~H,qA)- 

Jq 

(20) 

Again, as was the case for the LPA, the approx- 
imation at s = 2 provides information on all the 
n-point functions. This time, the n-point func- 
tions depend on a single momentum. They may 



be obtained as derivatives of the 2-point function, 
according to 

ri") (p, -p, 0, • • • , 0, 4>) = a^-^rf (p, 0), (21) 

which may be viewed as a generalization of 
Eq. Thus, for instance, the momentum de- 

pendence that remains within the 3 and 4-point 
vertices in Eq. (|19p is indeed that of the 2-point 
function itself. 

At this point an important subtlety appears, 
coming from the fact that the flow of the poten- 
tial (or of its second derivative) can be calculated 
either from Eq.®, in which Gk{q, (j)) is obtained 
from (HH) and ((TUl), or directly from Eq.([Tn]) at 
p = 0, since r^^^(0, (j)) = d'^Vk- If no approxima- 
tions were made, both results would be identical. 
However, once approximations are done, as it is 
the case here, both results do not coincide. 

At any order s of the BMW approximation 
scheme, the same ambiguity takes place for any 
correlation function P^"^ up to n = s — 1. Given 
the fact that the approximation is imposed only 

(s) 

on the flow equation of P^ and not on those 

(n) 

of the P^ with n < s, it is natural to com- 
pute these functions from their own flow equation 
(which is exact) and not from the flow equation 

of P^^'' (which is approximate). 

(s) 

One then subtracts from Tf. the parts of it 
that can be expressed in terms of lower order 
correlation functions, and perform the BMW ap- 
proximation in the equation for the difference. 
For s = 2, this amounts to computing the poten- 
tial from Eq.® and to implementing the BMW 

(2) (2) 

approximation on P^ '{p, (f>) — Tl. ' {0 , (f>) . 

The rationale behind this choice is that com- 
puting the flow of 9|Vfc from the equation for 

(2) 

rJc (P: 0) at p = would imply two approxima- 
tions: the equation for F^ (p, 0) is itself 
approximated, and propagators and vertices in 
its r.h.s. are also approximated. On the con- 
trary, Eq.® for the potential is formally exact 
and only the propagator used in it is approxi- 
mated. This general consideration can be made 
more concrete in the perturbative regime: The 
function F^ {p, (j)) obtained from Eq. (IT^ is one- 
loop exact and so is f[,^^ {p = 0, </>) = d'^Vk- When 
the corresponding propagator, computed from 
Eq. pH]) . is inserted in Eq.®, the obtained po- 
tential becomes two-loop exact. By generalizing 
the above subtraction procedure at higher orders 
similar perturbative considerations can be made: 
at order s = 2s' of the BMW scheme, the poten- 
tial computed from Eqs. (|9ll0p is (s' + 1)— loop 

exact, Tj, computed from Eq. pTI) is s'— loop ex- 
act, and so on. We thus expect that implement- 
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ing the BMW approximation only on the part 

(s) . . . 

of which is genuinely of order s, will have a 
decreasing impact on the lower order correlation 
functions as s grows. 



In practice, for s = 2, we rewrite 



rf ^(p, ^) = p' + Afe(p, ^) + (22) 

where Vk {(p) is obtained by solving ^ , and more- 
over, for numerical convenience (see below), the 
bare term has been extracted. The BMW ap- 
proximation is implemented only on the flow for 
lS.k ■ The equation for (p, (j)) can be deduced 
from (fT9| by subtracting its p = form: 



dt/^kip. P) = ^pMp, p) [uk{p) + Aiip, p)f - 2ph{p) uUp) ~ -h{p) [Kip. P) + 2pAi:(p, p)] , (23) 



with 



1 

T 

mlip) 
Uk{p) = 



(0,P) 



dpmlip) 



(24) 

(25) 
(26) 



and the symbol ' denotes the derivative with re- 
spect to p. 

In closing this section, let us mention that the 
relationships between the BMW scheme at order 
■3 = 2 and, on one the hand the large N expansion 
and, on the other hand the DE, are discussed 
respectively in Sections IV CI and IVIII 



IV. IMPLEMENTATION AT 
CRITICALITY 

In order to treat efficiently the low momentum 
region at criticality and, in particular, to cap- 
ture accurately the fixed point structure, we first 
introduce dimensionless and renormalized vari- 
ables, to be denoted with a tilda. We thus in- 
troduce a renormalization factor Z^, which re- 
flects the finite change of normalization of the 
field between the ultraviolet scale A and the scale 
k. Within the DE at O(V^) this factor describes 
the overall variation with k of the function Zk{(j)) 
in Eq.([n]). We define here Zk by 



Zi. 



dp^ 



(27) 



p=Po;p=po 



where pq and po ^^re a priori arbitrary. From Z^ 
we define the running anomalous dimension r\k 

by 



-kdk \nZk. 



(28) 



r 



Momenta are naturally rescaled according to 
p = kp. Other quantities are made dimension- 
less by dividing them by appropriate powers of k 
(and possibly conveniently extracting numerical 
factors). Thus we define 



p = Kdk^-^Z-'p, mUp) = Zk^ fhlip), 



-1 



Uk{p) = Zlk^-'K~'uk{p), 



1 „~, 



(29) 



where Kd is a constant originating from angular 
integrals. 



KT' =. 2 



We also set 

Gk{p,p) 
Jn{p,p) 



1 



d7r'^^^T{d/2). 



;Gk{p,p), 



Zkk^ 

Kd Jn{p,p)- 



(30) 

(31) 
(32) 



Instead of Rkil) it is convenient to work with 
a dimensionless cut-off function considered as a 
function of y = q^/k"^: 



Rkjq) 

q^Zk ' 



(33) 



Now, we note that as p — >■ at fixed fc, 
Ak{p,p) oc p^. This p^-dependence may gen- 
erate numerical instabilities in the equation for 
A (once transformed to dimensionless variables) . 
In order to avoid these, we found it convenient 
to introduce the renormalized and dimensionless 
2-point function Yk {p, p) : 



f 



Afe(p, p) 



Zfe(l + ffe(p,p)), (34) 



The function Yk{p, p) is a slowly varying function 
of p and its fiow equation is regular. This equa- 
tion is easily obtained from the flow equation for 
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A{p,p), Eq. It reads 

dtYk mi'^ + Yk) +pdpYk - (2 - d 
+2pp-^[{p^Y^ + Ukfj3- 
-hin/2 + pYl') . 



Vk)pYi 



ull: 



(35) 



The normalization condition ()27p is now ex- 
pressed as: 



^fc(p = Po,P = Po) = 0, Vfc. 



(36) 



The equation for (p, p) needs to be com- 
pleted by the flow equation for the dimensionless 
effective potential Vk{p) = k~'^Vk{p), or rather, 
the equation for its derivative Wk{p) = V(.{p), 
which is more convenient since Vk contains a triv- 
ial constant part which induces a numerical di- 
vergence. The equation for reads: 



dtWk{.p) = -i2-rjk)Wkip) + id-2+Tjk)pWl,ip) 

(37) 



Finally rjk in Eq. psp is implicitly determined 
by inserting the renormalization condition p6|) 
in the flow equation of Yfe(p, p) and evaluating 
the r.h.s. at p = po a-nd p = pQ. 

In principle, and if no approximations were 
performed, no physical quantity would depend on 
the choice of po and po j nor on the relatively free 
choice of the cut-off function r{y). In practice, 
the BMW scheme, as any approximation scheme, 
introduces spurious dependence on these choices. 
Below, we use the one-parameter family of cut-off 
functions: 



riy) = 



a 



cv -1 ' 



(38) 



and study the dependence of our results on po, 
Po, and a [l^. 

We numerically solved the flow equations from 
given initial conditions: 



WAip) -r/A2 
Ya{p,p) -0, 



{u/3)A'i-^Kdp 



(39) 



searching for the critical point by dichotomy on 
the initial parameters. We set u = uA'^~'^Kd- 

The numerical resolution is done on a fixed, 
regular, (p, p) grid, with < p < Pmax and 
< p < Pmax- With our choice of cut-off func- 
tion, the contribution of the momentum interval 
q £ [4, oo] to the integrals /„ and J3 is extremely 
small and thus we neglect it by restricting the in- 
tegration domain to g G [0,4]. When computing 
the double integrals Jsip, p), we need to evaluate 
Y for momenta p + q beyond Pmax- In such cases, 
we set Y{p > Pmax) = ^'(Pmax), an approxima- 
tion checked to be excellent for p,„ax > 5. To 



0.02 
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FIG. 1. (color online) Example of the dichotomy pro- 
cedure used for reaching the fixed point (A'^ = 1, 
d = 3, u = 6 X 10"^po = 0,po = 0, a = 2.25). 
The plot shows the running anomalous dimension rjk 
as a function of t = ln(fc/A). Each curve corresponds 
to a different initial value of r (see Eq. (|39p ). The red 
dashed line indicates the estimated asymptotic value 
77 ~ 0.03943. 



access the full momentum dependence, we also 

(2) 

calculate r)^, {p, p) at a set of fixed, freely cho- 
sen, external p values. For a given such p, p/k 
is within the grid at the beginning of the fiow. 
This is no longer so when k < p/pma.x] then, we 
switch to the dimensionful version of Eq. psp . 
and also set J3(p, p) = G(p, p) J2(0, p), an excel- 
lent approximation when p > kp^^^^. 

We found that the simplest time-stepping (ex- 
plicit Euler), a finite-difference evaluation of 
derivatives on a regular (p, p) grid, and the use 
of Simpson's rule to calculate integrals, are suf- 
ficient to produce stable and fast-converging re- 
sults. For all the quantities calculated, the con- 
vergence to at least three significant digits is 
reached with a (p, p) grid of 50 x 60 points and 
elementary steps 5p = 0.1 and 5p = 0.1. With 
such a grid, a typical run takes a few minutes 
on a current personal computer. The step in 
t = logfc/A is Ai = 10"'', and the flow is run 
down to i ~ —20. In order to find the fixed point, 
we performed a simple dichotomy procedure on 
the bare mass m\ = r/A^ at fixed u, by studying 
the flow of Wk (0) . Fig. [T] iUustrates the flow of 
Tjk as one approaches the fixed point. 



V. RESULTS AT CRITICALITY 

Although the main goal of the BMW method is 
to provide access to full-momentum dependence, 
it can of course also be used to compute criti- 
cal exponents and other zero-momentum quanti- 
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FIG. 2. (color online) Running exponent {d = 3, 
N = 1, a = 2.25, po = and po = 0). Each curve cor- 
responds to a different initial value of r (see Eq.(|39}). 
Inset: the exponential approach to the asymptotic 
exponent is used to estimate uj ~ 0.784. 



ties [1^. In this section, where we return to the 
0{N) models (with general N), we provide de- 
tails on the calculation of the critical exponents 
and check their robustness with respect to vari- 
ations of the different parameters of the method 
such as the numerical resolution, the choice of 
the cut-off function and the location of the nor- 
malization point (pojPo). 

Since we focus here on the regime of small mo- 
menta, it is convenient to take as initial condi- 
tion a value of the dimensionless coupling u not 
too small compared to 1 in order to initialize the 
flow far from the Gaussian fixed point u = and 
thus to approach quickly the infrared fixed point. 
This is useful, not only because of the short- 
ened time needed to reach the critical regime, but 
also because otherwise, due to the 16-digit pre- 
cision used, our dichotomy procedure does not 
allow for an accurate determination of the fixed 
point directly from initial parameters. The re- 
sults to be presented below have been calculated 
for u = 6 X 10-^ /N. 



A. Numerical extraction of critical 
exponents 

The anomalous dimension rj comes out of the 
solution of the flow equations, which provide a 
direct estimate of rjk (Fig. [1]). It can also be 
extracted from r[,^p(p, p = 0) oc p'^^^ at small 
momentum with the same result, although this 
is a much less practical way. 

In the vicinity of the fixed point, the behavior 
of any dimensionless and renormalized quantity, 
such as the dimensionless mass, is as follows (re- 



call that < = Infc/A < 0) 

rnl = fhl+ m\cr^^ + 77126"^* -|- msc""* -f . .(40) 

with the universal critical exponent v describing 
the departure from the critical surface, and the 
correction to scaling exponents w, W2, ■ • • describ- 
ing the initial approach to the fixed point. 

In practice, we use the flow of the mass (the 
flow of ?7fe could also be used) to extract v and w 
[20| . We explore successively regions of t- values 
where one of the exponentials in the equation 
above dominates. For instance, for t negative 
enough, we write: 

log 19(771^1 '-^ [-constant (41) 

to find V. To extract w, we choose |i| large enough 
but not so large as to leave the vicinity of the 
fixed point. We then write 

log |9tm^| ^ Lai + constant. (42) 

Notice that away from the fixed point, the expo- 
nents thus determined depend themselves weakly 
on t since, strictly speaking, pO)) holds only in 
the infinitesimal vicinity of the fixed point. We 
thus obtain only (slowly) running exponents i/fc 
and ujk. In practice, these exponents are calcu- 
lated by taking the i-derivative of Eqs. (|41l42p . 
The procedure is then repeated for a set of ini- 
tial conditions that bring the system closer and 
closer to the critical point. The estimates of i/fe 
and Wfe saturate to their fixed point values re- 
flected in the plateau seen in Fig. [2] for (w is 
the most difficult exponent to determine numeri- 
cally). Given such curves, one can further extract 
even more accurate estimates from the (exponen- 
tial) approach to the asymptotic plateau values, 
see inset of Fig. [21 

With this method we could, in principle, ex- 
tract exponents with almost arbitrary numerical 
accuracy. In practice, however, only a few dig- 
its are significant: our results suffer indeed from 
an uncertainty related to the choice of the cut- 
off function (see next subsection); besides, it is 
not necessary to present results with an accu- 
racy that far exceeds the deviation from those 
with which they are compared. 

B. Dependence on renormalization point 
and regulator 

Although as explained above the values of the 
critical exponents should in principle depend nei- 
ther on the normalization point (poiPo) nor on 
the shape of the cut-off function Rk{q) this is 
no longer the case once approximations are per- 
formed. 
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In practice, we apply the "principle of mini- 
mal sensitivity" , searching for a local extremum 
of the physical quantities under study [l^ [llj in 
a "reasonable" subspace of values taken by a (the 
parameter of our cut-off function (pSl) ). po, and 
po- It is then expected that the corresponding 
values are "optimal" in the sense that they show, 
locally, the weakest dependence on the above pa- 
rameters. 

Here, we first notice that at fixed a and po, 
the dependence of our estimates on pQ is much 
weaker than that found by varying a and po. 
Fig 131 shows the variation of the anomalous di- 
mension r/ with a for two typical values of N in 
three dimensions. As in all other cases studied, 
we observe the existence of a unique extremum. 
In the following, we always use these extremum 
values to report our best estimates for the crit- 
ical exponents. Note that we do not show the 
variations of the exponents with po as they can 
be shown to be equivalent to those with a [22j . 




a 



FIG. 3. (color online) r; as a function of the cut- 
off parameter a, for A'' = 1 and = 10 (d = 3, 
Po = 0,po = 0). 



quoted remain stable when a varies in the range 
[a* -1/2, a* + 1 /2] . The quality of these numbers 
is obvious: our results for u agree with previous 
estimates to within less than a percent, for all N; 
as for the values of 77 and w, they are typically 
at the same distance from the Monte-Carlo and 
high temperature series estimates (for instance, 
hr N = l,u = 0.6298(3) [H) as the results from 
resummed perturbative calculations. Our num- 
bers also compare favorably with those obtained 
at order in the DE scheme [l9l |. 



In the limit of large N, the BMW scheme be- 
comes exact for the 2-point function for s > 2 
0, [1]. This generalizes the fact, shown in p^ . 
that the LPA (s = 0) is exact in the large N 
limit for the effective potential. It can be veri- 
fied from the tables 1 to 4 that the large N limit 
values T] — 0, V = 1 and u ~ 1 are approached 
for large values of N. 

We can also perform a 1 /N expansion [2^, [2^ . 
This was already done in where the BMW 
scheme was further approximated by the use of 
LPA propagators. As the LPA becomes exact in 
the large N limit, these results are unchanged at 
first order in 1/iV, except for the use of another 
type of regulator profile. An analytical study of 
the BMW equations in this limit provides the fol- 
lowing values for the critical exponents at order 
1/iV: ry = 0.23/7V and i/ = 1 - IMi/N, to be 
compared with the exact results [1^ r] = 0.21 /N 
and I' = 1 - 1.081 /Af. In [ll| the use of another 
regulator profile allowed us to obtain somewhat 
better results for 77 in this limit: 77 = 0.25/A^. 
Notice that all these analytical results are recov- 
ered in our numerical solution for large values of 
N (notice in fact that terms of order 1/A^^ are 
very small already for > 4). 



C. Results for the critical exponents 

We now present our results for the critical ex- 
ponents of the scalar 0{N) models in d = 3. 
They have been obtained with a two-dimensional 
grid in p and q with rip = 51 points in the p- 
direction, iiq — 60 points in q and with qmax — 4, 
Pmax = 6, and a pmax = 5N. Tables HI HIl and Hill 
contain our results for the critical exponents 77, 
v and w, together with some of the best esti- 
mates available in the literature, obtained either 
from Monte Carlo or resummed perturbative cal- 
culations (that we refer to as field theory (FT)). 
Our numbers are all given for the optimal val- 
ues a* of the cut-off parameter, and the digits 



The two-dimensional case, for which exact re- 
sults exist, provides a very stringent test of the 
BMW scheme. We focus here on the Ising model 
A^ = 1 which exhibits a standard critical behav- 
ior in d = 2, and the corresponding critical expo- 
nents. Notice that the perturbative method that 
works well in d = 3 fails here: for instance, the 
fixed-dimension expansion that provides the best 
results in d = 3 yields, in d = 2 and at five loops, 
77 = 0.145(14) 27] in contradiction with the ex- 
act value 77 = l/4[? ]. We find instead 77 — 0.254, 
V — 1.00 in excellent agreement with the exact 
values r/ = l/4, z/ = l. A more detailed study of 
0{N) models in d = 2, at and out of criticality, 
will be presented in a separate work. 
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TABLE I. Results for the anomalous dimension in d = 3, compared with results obtained within the DE at 
order O(V^), field theory (FT) and Monte Carlo (MC) methods 



iV 


BMW 


DE 


FT 


MC 





0.034 


0.039[28l 


0.0272(3) [29] 


0.0303(3)) [30] 


1 


0.039 


0.0443 [19] 


0.0318(3) [29] 


0.03627(10) [31J 


2 


0.041 


0.049 [28] 


0.0334(2) [29] 


0.0381(2)[^ 


3 


0.040 


0.049[28] 


0.0333(3) [29] 


0.0375(5) [33] 


4 


0.038 


0.047[28] 


0.0350(45) [34] 


0.0365(10) [35] 


10 


0.022 


0.028[28] 


0.024 [36] 




100 


0.0023 


0.0030[28] 


0.0027 [26] 




0{1/N) 


0.23/iV 


0.27Q/N [26] 





TABLE II. Results for the critical exponent v in d = compared with results obtained within the DE at 
order O(V^), field theory (FT) and Monte Carlo (MC) methods 



iV 


BMW 


DE 


FT 


MC 





0.589 


0.590[28] 


0.5886(3) [29] 


0.5872(5) [37] 


1 


0.632 


0.6307[19] 


0.6306(5) [29] 


0.63002(10) [31J 


2 


0.674 


0.666 [28] 


0.6700(6) [29] 


0.6717(1) [32] 


3 


0.715 


0.704 [28] 


0.7060(7) [29] 


0.7112(5) [33] 


4 


0.754 


0.739 [28] 


0.741(6)[34] 


0.749(2) [35] 


10 


0.889 


0.881[28J 


0.859 [36] 




100 


0.990 


0.990 [28] 


0.989[26] 




0{1/N) 


1 - l.QM/N 




1 - 1.081/iV [26] 





TABLE III. Results for the correction to scaling ex- 
ponent uj in d = 2) compared with results obtained 
within the BMW method, field theory (FT) and 
Monte Carlo (MC) results 



iV 


BMW 


FT 


MC 





0.83 


0.794(6) [29] 


0.88 [30] 


1 


0.78 


0.788(3) [29] 


0.832(6) [31J 


2 


0.75 


0.780(10) [29] 


0.785(20) [32] 


3 


0.73 


0.780(20) [29] 


0.773 [33] 


4 


0.72 


0.774(20) [34] 


0.765 [35] 


10 


0.80 






100 


1.00 







D. The function F'^' at criticality and 
further tests at intermediate and large 
momenta 

We now study the momentum dependence of 
the two-point function at criticality. In dimen- 
sion three, the bare couphng constant u has 
the dimension of a momentum and thus sets 
a scale (the Ginzburg length: ~ u^^). 

There are typically three momentum domains for 
r(2)(p,p = 0) PO: 

(i) the infrared domain defined by p <C m where 
r(2)(p) ^ -uV"''- We show in Fig. Hthat this 
behavior is well reproduced by our solution of 
the flow equation. To clearly see this regime on 
a large range of momentum we have integrated 
the flow with a bare value of u not too far from 




FIG. 4. (color online) The ratio of the 2-point func- 
tion r(2' (p, 0) and of ^ at criticality as a func- 
tion of p/A (d = 3,Ar = 2,a = 2,po = 0,po = 0). 
The normalization has been chosen so that this ratio 
starts close to 1 at small p. The bare dimensionless 
coupling is S = 6.10^^/A'^. 



the value of A: u = 6.10"V^- 

(ii) the ultra-violet domain defined hy p ^ u 
(and h.^ p) where T^'^'>{p) can be studied per- 
turbatively and is found to behave at two loops 

as r(2)(p) - p2 (CAr/967r2)u2 log£/u (with 

= {N + 2)/3). It was shown in [iJl that in 
the BMW approximation, and at large momenta 



10 



10 



10" 



(2-11) 

P 

(l/96n) u log( p/u) 



have calculated the quantity 




10 10 



10"' /A 10° 
p/A 



at 



FIG. 5. (color online) The difference T'^^ (p, 0) 
criticality as a function of p/A, compared with its ex- 
pected ultra-violet behavior ~ logp/w (d = 3, A*' = 
2,a = 2,po = 0, po ~ 0). The infrared p^"'' behav- 
ior is also shown (see text). The bare dimensionless 
coupling is ti = 10~^/A'^. 



TABLE IV. Results for the quantity c defined in the 



text. 


A' 


BMW 


lattice 


7 loops [42] 


1 


1.15 


1.09(9)[43] 


1.07(10) 


2 


1.37 


1.32(2) [44| 


1.27(10) 






1.29(5)[441 




3 


1.50 




1.43(11) 


4 


1.63 


1.60(10)[43J 


1.54(11) 


10 


2.02 






100 


2.36 







A(p, p = 0) behaves as In (p/u), more precisely, 
dAk=o{p,0) 



d\p\ 



C 



N 



2|p| 



dtRk{l)Gtil)Gt,{q). 



(43) 

The u^logp/u behavior is thus retrieved, see 
FiglSl with a prefactor that however depends on 
Rk{q)- With the exponential cut-off function, 
Eq. (l38t , the prefactor can only be calculated nu- 
merically. We have studied its dependence on 
a and shown that there is an extremum around 
a ~ 5 where the difference with the exact re- 
sult is about 8%. Of course, this UV behavior 
shows up only if the bare coupling u is suffi- 
ciently small compared to A. We have chosen 
u = 10~^/N to have a large UV domain where 
this behavior is clearly seen. Note that at small 
p, r^^^ (p) — p^ p^~'^ which is visible on FiglS] 
although this regime is approached very slowly. 

(iii) the cross-over between the infrared and 
ultra-violet domains. This regime of momentum 
is visible on both figures 2] and [5] for p ~ u. 

For purposes of probing the intermediate mo- 
mentum region between the IR and the UV, we 



d^p 



r(2)(p) 



(44) 



which is very sensitive to the cross-over regime: 
the integrand in Eq. (|44| is peaked at p ~ 
(Nu)/ 10 [38;]. For this reason, the calculation 
of c has been used as a benchmark for non- 
perturbative approximations in the 0{N) model. 

In the 0{2) case and for d — 3, this quantity 
determines the shift of the critical temperature 
of the weakly repulsive Bose gas [111 (notice that 
c is not defined for d = 2). It has thus been much 
studied recently using various methods, even for 
other values of N. In particular, the large N 
limit for this quantity has been calculated an- 
alytically and found to be c = 2.3 [i^. In this 
work, we have found the values for c for some rep- 
resentative values of N. Our results, compared 
to the best ones available in the literature (with 
their corresponding errors when available), are 
presented in Table IIVI For all values of N where 
lattice and/or 7-loops resummed calculations ex- 
ist, our results are within the error bars of those 
calculations (and comparable to those obtained 
from an approximation specifically designed for 
this quantity 0, |4l|), except for = 2, where 
very precise lattice results are available. In the 
large N limit, one can see that our results dif- 
fers from the exact value by less than 3%. Notice 
that the large N behavior of the quantity c is in 
fact of order 1/A^ [i^, which as we have seen is 
not calculated exactly at this level of the BMW 
approximation. 

Altogether, we can see that the BMW method 
is able to reproduce the correct behavior of the 
2-point function at criticality in all momentum 
regimes. Note in particular that this is not the 
case of conformal field theoretical methods that 
are only able to capture at criticality the confor- 
mally invariant p^~^ behavior but that can re- 
produce neither the ultra-violet behavior, corre- 
sponding to u <^ p <^ A, nor the cross-over be- 
tween the infrared and ultra-violet regions, cor- 
responding to p ~ u. 



VI. SCALING FUNCTIONS 

As an approximation of the NPRG, the BMW 
scheme allows us to investigate all momentum, 
temperature and external magnetic field regimes, 
and is not restricted to the long distance physics 
at criticality. A particularly interesting, and 
a priori difficult regime is the critical domain, 
where the correlation length is large but finite. 
In this case, an appropriately rescaled two-point 
function shows a universal behavior. As the 
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BMW approximation allows for the calculation 
of genuine momentum-dependent quantities, the 
calculation of this scaling function and its com- 
parison with the best available theoretical results 
from the literature and with experimental data 
represent one of the most stringent tests of the 
approximation. 

In this work, we consider the case N = 1 rel- 
evant, for instance, to describe the critical be- 
havior of fluids near the liquid-gas critical point. 
Near this point and for p <^ ~ m one expects 
the general scaling behavior 



(45) 



with, by definition, G^^-' the density-density cor- 
relation function, — r'^^^(p = 0) the com- 
pressibility and — k'^m'l, with fc — 0, the 
correlation length that diverges close to critical- 
ity with the v critical exponent. Here ± refers 
to the two phases, above and below the critical 
temperature respectively. The functions g±{x), 
normalized so that 



(46) 



are universal. Their limiting behavior is well 
known. For small x they are well described by the 
Ornstein-Zernicke (mean-field) approximation; 



goz{x) 



1 



1 + x^ 



(47) 



The corrections to the Ornstein-Zernicke behav- 
ior are usually parameterized as [46j 



n=2 



(48) 



The above behavior of g±{x)~^ is a priori valid 
only for x < 1 but since the coefficients c„ 
are very small, it turns out that the Ornstein- 
Zernicke approximation is actually valid over a 
wide range of x values, as we shall see later. For 
large x (that is, f ^ p~^) the scaling functions 
show critical behavior with an anomalous power 
law decay 



9±{x) 



(49) 



which allows for the experimental determination 
of the exponent 77. This expression also allows for 
corrections, as given by Fischer and Langer 



(50) 



Different approximate results for the universal 
scaling functions exist in the literature, obtained 



either by Monte Carlo methods '46'] , or by the use 
of an analytical ansatz, interpolating between the 
two know limiting regimes (|48p and (|49p , using e 
expansion results (the Bray approximation [48|). 
Experimental results from neutron scattering in 
CO2 near the critical point also exist [49j. 

In Bray's interpolation for the high tempera- 
ture phase one assumes g^(x) to be well defined 
in the complex x^ plane, with a branch cut in 
the negative real x^ axis, starting at x^ — ~t\^ 
where r\ = dM^^^^^ = 9Sm, following the theo- 
retical expectation that the singularity of g+{x) 
nearest to the origin is the three-particle cut 
[isl . [50l |. The parameter Mgap is the mass gap 
of the Minkowskian version of the model. For 
the 0"* theory, it is known that the difference be- 
tween the mass gap and is very small and re- 
placing one by the other corresponds to an error 
which is beyond the accuracy of our calculation 
[46| . Then Bray's ansatz in the high temperature 
phase (the only phase studied in the following) 
reads: 

_i 2 sin7r?7/2 



X I du (u) 

lr+ 



Sm 



Sm + x"^ 



(51) 



where F+(u) is the spectral function, which sat- 
isfies F+(-|-oo) = 1, F+(w) = for u < r_|., and 
> for M > r+. On top of this, one must 
impose g~^{Q) = 1, which fixes the value for Ci- 
One must then specify F+{u). Bray pro- 
posed the use of a spectral function with the ex- 
act Fischer-Langer asymptotic behavior, of the 
type 



P+{u) — Q+{u) cot ^nrj 

p+{uY+Q+{ur 



(52) 



where 



^+(«)l + ^cos^ 
^ ,^C+ . ttC 



— TT + —rr cos — 



C+ TT 

— ri— sm — 



(53) 



with C = (1 — a) /v. This definition contains a 
certain number of parameters. On top of the 
critical exponents, which can be injected using 
either the BMW values or the best available re- 
sults in the literature, one must also fix 5^^, C2 



and Cq 



For S^.j one can use the best estimate 



in the literature, given by the high temperature 
expansion of improved models [51|. Bray pro- 
posed to fix C2 + to its e-expansion value 
C2 + — —0.9, and to then determine by 
requiring F+^b{'I'+) = 0- These conditons allows 
for a little parameter tuning, by adjusting the 



12 




FIG. 6. (color online) The 2-point scaling function 
g{p£,) as a function of a; = in the high-temperature 
phase {d = S, N = 1). Solid blue line: BMW re- 
sult. Red squares: Experimental results of [i^. In- 
set: same data with logarithmic scales. 



relative weight of the and parameters. 
When comparing our results with Bray's ansatz, 
we shall use this freedom. We now turn to the 
scaling function computed by the BMW method. 

In terms of the variables used in this paper, we 
find that 



(54) 



when fc — > 0. In this work, for purposes of com- 
parison with existing results, we have only com- 
puted the high temperature scaling function. We 
have performed the calculation for different val- 
ues of the correlation length (and hence of the 
reduced temperature). When plotted, one can 
indeed see perfect data collapse for different val- 
ues of which is the first non trivial test of the 
quality of our results for the scaling function. 

In Figini we plot the BMW scaling function 
together with the experimental results from ref- 
erence [49|. Due to the small values taken by 
the coefficients c„ and the critical exponent rj in 
d = 3, the Ornstein-Zernicke behavior dominates 
even beyond = 1. In order to measure the 
deviation from this behavior, one usually makes 
use of the auxiliary function 



h{x) = log 



goz{x) 



(55) 



In Figl7]we plot this function together with the 
experimental results from [4§| and the results 
from the Bray ansatz for two "extreme" choices 
of the C2 and parameters. One can there 
see that the BMW approximated result compares 
very well with all these results. In particular, it 



is in between the results obtained from the two 
Bray ansatz considered. 

Let us mention that even with large system 
sizes, the Monte Carlo results suffer from signif- 
icant systematic errors for larger than typi- 
cally 5 to 10. This probably comes from the fact 
that the universal behavior of the structure fac- 
tor shows up only when ^ and the separation I 
between the spins at which we calculate the cor- 
relation function are large compared to the lat- 
tice spacing and small compared to the lattice 
size: even for lattice sizes of a few hundreds of 
lattice spacings this leaves only a small window 
of useful values of ^/Z [46| . 



h(x) 
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FIG. 7. (color online) Deviation of the scaling func- 
tion to its trivial Ornstein-Zernicke form, Ea. (|55|) . 
The dotted and dashed lines correspond to two "ex- 
treme" choices of the parameters Ci and C3 of 
Bray's ansatz. Dotted line: Ci = 0.924, C2 = 1.8, 
C3 = -2.28. Dashed hne: Ci = 0.918, C2 = 2.55, 
C3 = -3.45. 

On top of these results we can also compare 
results for the values of the coefhcients and 

Ct- The resuhs for BMW are c+ 4.5 x 10"" 

to be compared with the IHT best estimate [l^l 
4 = -3.90(6) X 10-^ whereas for C+ BMW 
yields — 0.914, to be compared with the e- 
expansion result — 0.92. 

We conclude this section by noting that (i) the 
structure factor encompasses much more infor- 
mations on the universal behavior of a model 
than the (leading) critical exponents (that are 
moreover difficult to measure experimentally), 
(ii) Bray's ansatz, although powerful, depends on 
two parameters C2 and C3 that are poorly deter- 
mined perturbatively as well as on two critical 
exponents, (iii) the present state of the art of the 
Monte Carlo simulations is by far insufficient to 
compute reliably the structure factor in the inter- 
esting region of momentum where is large, (iv) 
the BMW method leads to a determination of the 
structure factor that has no free parameter once 
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a choice of regulator has been made (possibly in- 
volving an optimization procedure as described 
in Section IV Bp . The results above, summarized 
in Figl?! suggest that the BMW method leads to 
an accurate determination of the structure factor 
in the whole momentum range while the experi- 
mental results seem to suffer at small momentum 
from systematic deviations. 



VII. RELATION WITH THE 
DERIVATIVE EXPANSION 

The validity of the DE is rarely questioned, 
satisfactory results being taken as an a posteriori 
check. We show now that the BMW approach 
allows for a deeper understanding of its range of 
applicability, and of some of its peculiar features. 

The ansatz defining the order of the DE (see, 
for instance, Eq. p^ for the order 2) is used to 

(i) define the quantities to be determined, 
which, in the case of order 2, are the effective 
potential and the field normalization, both func- 
tions of the (constant) field 

Vfe(0) = ^r,[0(x)]|^(^^^^ 

(56) 

^z,(^) = a,. (ri^) [;.,</>]) 

(ii) compute the n-point functions r^"-* and the 

propagator Gk — (F^^^ + Rk)^^ that enter the 
right hand sides of the flow equations of , Z^^ 
etc. 

In short, the DE projects the functional Tk[(f)] 
on a polynomial expansion in powers of the 
derivatives of the field, the expansion coefficients 
being field dependent. In Fourier space, the DE 
amounts to a polynomial expansion of the n- 
point functions F^,"'' (pi , • • • ^Pmfp) in powers of 
the momenta Pi, around vanishing momenta (see 
for instance Eq. dTSt V At this point, it is use- 
ful to introduce a distinction between external 
momenta, the momenta that appear in the n- 
point function T^j^\pi, ■ ■ ■ ,pn,(f>) whose fiow is 
being considered, and the internal momentum, 
denoted by q, appearing in the n-point functions 
in the r.h.s. of the corresponding flow equation 
and which is integrated over. In contrast to what 
is done in the BMW approximation, in the DE 
no distinction is made between these two sets of 
momenta, which can lead to inconsistencies. For 
instance, in the fiow equation for Zk{(f>) at order 
2, the product F^ '{p,q, -p-g)F^ '{-p, -q,p+q) 
(see Eq. pT]) ) leads to four terms of order four: 
{p^Y ,p'^q^ ,{p.qY ,{q^Y , that, in a strict expan- 
sion to this order, should be neglected (note that 



this is not what is usually done in the DE con- 
text). In fact, since Zk{4>) is already the coeffi- 
cient of the p^ term in the expansion of F^j^^ (p, 0), 

any dependence of F^'^^ (and of F^'''' ) on the inter- 
nal momentum q should be neglected in dkZk{(j)) 
at this order of the DE. Since the BMW approx- 
imation at order s — 2 precisely consists in set- 
ting q = in F^'^-' and F^*'' in the flow equation 

of F^ , we conclude that at this order the BMW 
approximation contains all terms of the DE at 
order V^. 

The BMW approximation, that disentangles 
the roles of the internal and external momenta, 
differs deeply from the DE precisely on the point 
explained above: as the DE, it takes advantage 
of the fact that the internal momentum is cut-off 
by dkRkiq) in order to expand in powers of q/k 
(in fact only the leading term, g = 0, is retained), 
but does not rely on the smallness of the external 
momenta. 

In fact, the natural expansion parameter of the 
DE is the ratio p/k or p/m, whichever is small- 
est, where m is the smallest of the masses that 
may appear in the problem considered: When k 
is much larger that all masses, these can be ig- 
nored and p/fc is the expansion parameter. When 
k becomes smaller that the smallest mass, the 
flow essentially stops and the expansion parame- 
ter becomes p/m in the limit fc — ^ 0. Thus, it is 
plausible that the DE performed as a power se- 
ries in p/k in a critical theory (m = 0) possesses 
a radius of convergence of the same order as the 
DE performed as a power series in p/m in a mas- 
sive theory at fc = 0. In this last case, the radius 
of convergence is known for = 1 in dimension 
three (' [isl . [50|): It is 3 in the symmetric phase 
and 2 in the broken phase [? ] . 

The above arguments suggest that the DE 
is not able to describe fc-dependent correlation 
functions with external momenta higher than 
typically 3max(fc,m). In particular, in the crit- 
ical case where massless modes are present, the 
DE is only suited for the calculation of physical 
(that is at fc = 0) correlation functions a.t p ~ 0: 

(2) 

The anomalous momentum behaviour F)j^q(p) ~ 
p^^^ , valid at small p, will not emerge at any or- 
der of the DE. Of course, this does not mean that 
the anomalous dimension cannot be determined 
within the DE, as one can exploit general scal- 
ing relations and the fact that the anomalous di- 
mension enters also quantities that are defined at 
zero momentum. Thus, for instance, "q can be es- 
timated from the fc-dependence of the normaliza- 
tion factor Zk ~ k~^ (or alternatively from the 
large field behavior of the fixed-point dimension- 
less effective potential). In contrast, the BMW 
approximation correctly captures the anomalous 
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scaling of T\,J^q{p) at small p, and this is a di- 
rect consequence of the fact that no expansion in 
external momenta is performed [? ]. 

The origin of the difficulties of the DE is that 
it does not have good decoupling properties in 
the momentum range p ^ k. The decoupling 
property, crucial for universality, means, on the 
example on the 2-point function, that (p) be- 
comes almost fc-independent when k <^ p and 
that therefore T^-^loip) ^ ^^klpip)- One could 
thus naively expect that external momenta {pi}, 
i — 1, . . . ,n play the role of infrared regulators in 
the flow of r[,"-'({pi}) and that when k < pi,\fi 

the flow of r^"-* (almost) stops. In fact, in flow 
equations, external momenta play, at best, the 
role of infrared regulators when all momenta in- 
volved (external and internal) are not in an ex- 
ceptional configuration. The problem is that 
even when the external momenta are not excep- 
tional, the integral over the internal momentum q 
in the flow equation of F^"^ involves vertex func- 
tions (r[,"^^-' or r^"^^-*) in exceptional configura- 
tions. Depending on the approximation scheme, 
this can spoil the decoupling property that, un- 
doubtly, should hold for the (physical, that is 
k <^ p) correlation functions themselves when 
they are evaluated in non-exceptional configura- 
tions. The difficulty is therefore to devise an ap- 
proximation scheme that satisfies the decoupling 
property. While this is the case of the BMW 
scheme it is neither of perturbation theory nor of 
the DE. One can nevertheless try to extract from 
the DE the gross behavior of F^^^ (p) (and of the 
other functions) by stopping by hand the flow 

(2) (2) 

at k = p and identifying r^^p(p) with F^^p (p) . 
This idea has been explored in [§| (see also p^)- 
The resulting correlation functions roughly show 
the expected momentum behavior, but as ana- 
lyzed in detail in [1] , it does not seem possible to 
extend this first qualitative analysis and to ob- 
tain quantitatively precise correlation functions 
without having recourse to BMW. 

To gain further insight into the validity of the 
DE, we may consider a simple analytical repre- 

(2) 

sentation of the function F^^g (p, (f>) determined 
with the BMW approximation at order s = 2, 
which, as we have shown, is very close to the 
exact 2-point function over the whole momen- 
tum range. The following formula (inspired by 
Eq.(2.33) of (3) 

Tf\p,p) = Ap^p^ + bk' + b'Miip))-'"'^' 

(57) 

where A, b, and b' are independent of p and p, 
and (p) is a function homogeneous to a square 



mass, provides a good fit of the BMW results 
when k, as well as M'^{p), are very small com- 
pared to the ultra-violet cut-off A. This formula 
encompasses the two different regimes that char- 
acterize the behavior of Tj, (p) at small p: First, 
for p small compared to A and large compared 
to k and to the mass, it yields F^^^(p) ^ p'^~^'' , 
with rjk the running anomalous dimension. Thus, 
the critical behavior is captured for k sufficiently 
small for 77^ to be quasi-stationary and (almost) 

equal to 77. Second, for p small compared to ei- 

(2) 

ther Mfe or fc, one can expand F^, (p, 0) in powers 
of p2/(fc2 -I- M|) and get: 

rf)(p,p)=A' (fc2 + 6"Af|(p))-"^/^ 

Xp^l + hMp)p' + f2Ap)p' + ---) 

+ Vl + 2pVi:. (58) 

This is the kind of ansatz considered by the DE 
and it illustrates how the anomalous dimension 
can be extracted from the fc-dependence of the 
coefficient of the p^ term in the running action 
i- 

Finally, let us stress that the above remarks, 
while they provide some justification for the DE 
and in particular specify the conditions for its 
validity, are not sufficient to prove convergence, 
which may be strongly affected by the regulator. 
In particular, one may expect systematic errors 
in cases where the range of the cut-off function 
Rk{q) is not smaller than the natural radius of 
convergence of the DE. Notice however that at 
least for TV = 1 in d = 3, the smallness of the c„ 
coefficients in Eq. suggest that even at low 
order the DE should be able to capture the low 
momentum physics. An in-depth study of this 
issue will be presented in [22j . 



VIII. CONCLUSIONS 

In this paper we have presented the complete 
numerical implementation of the BMW approx- 
imation scheme that allows for a solution of the 
NPRG flow equations keeping the full momen- 
tum dependence of the 2-point function. At the 
level considered in this paper, this amounts to 
solve two coupled equations for the effective po- 
tential and the 2-point function. These equa- 
tions can be solved by elementary numerical tech- 
niques. 

We have considered applications to the 0{N) 
models, mostly in dimension d — 3. An accurate 
momentum dependence of the 2-point function 
has been obtained from the low momentum crit- 
ical region to the high momentum, perturbative. 
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region (such a region exists when the dimension- 
ful bare couphng is small compared to the ultra- 
violet cutoff). In particular, the critical expo- 
nents are accurately determined as was already 
reported in jl2t] . The additional results presented 
in this paper concerns the scaling functions which 
probe a different aspect of the momentum depen- 
dence of the 2-point function in the vicinity of the 
critical point. We have considered more specif- 
ically the scaling function for the case iV = 1 
above the critical point and have shown that it 
is in excellent agreement with the best available 
theoretical estimates. Interestingly, these esti- 
mates, including ours, differ significantly from 
the experimental data at small momenta. These 
scaling functions, which are difficult to obtain 
with other more conventional techniques, includ- 
ing Monte Carlo simulations, come out directly 
from the 2-point function obtained by solving the 
flow equations. 

Another information of physical interest which 
is also contained in the 2-point function that we 
compute is its field dependence. Thus a natu- 
ral application of the present method could be 
the investigation of the 0{N) models in the pres- 
ence of an external magnetic field. We could also 
contemplate extracting from the 2-point function 
information about possible bound states [5^ . Fi- 
nally, we note that the BMW method paves the 
way towards understanding a variety of situa- 
tions where the momentum structure plays a cru- 
cial role. For instance, a method similar in spirit 
has been applied successfully to the determina- 
tion of the fixed point structure of the Kardar- 
Parisi- Zhang equation [13, [11] and to the calcu- 
lation of the spectral function in a Bose gas [s^l . 
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Appendix A: Notation and conventions 

By taking successive functional derivatives of 
Tk[(j)] with respect to (pix), and then letting the 
field be constant, one gets the n-point functions 



S(l){xi) . ..54>{Xn) 



(Al) 

il>{x)=ct> 

in a constant background field 4>. Since the back- 
ground is constant, these functions are invariant 
under translations of the coordinates, and it is 
convenient to factor out of the definition of their 
Fourier transform the (5-function that expresses 
the conservation of the total momentum. Thus, 



with the usual abuse of notation, we define the 
n-point functions t')^\pi, . . . ,pn,4>) as: 

{2^f ri")(pi,...,p„,0)^ 



d^xi...d''xne'^^^^^'Tf\x^,--- ,x„,0). 

We use here the convention of incoming 
momenta, and it is understood that in 



,Pm4>) the sum of all momenta van- 



(n) 

ishes, so that F^ is actually a function of n — 1 
momentum variables (and of (j)). Notice that we 
use brackets for functional, e.g. Ffe[^], and paren- 
thesis for functions, e.g. F^"''(pi, • • • ,p„, 0) when 
4> is uniform. For the 2-point function evalu- 
ated in a uniform field configuration, which effec- 
tively depends on a single momentum p, we often 

(2) 

use the simplified notation Fj^, {p, (f) in place of 
Ff (p, 



Appendix B: Extension of BMW 

In the approximation BMW with s = 2, 
we make the following substitutions in the 

(2) 

r.h.s. of the flow equation for Fj, (p): 



> F^^^(p,-p,0,0), and 
r), ' {p,0, —p), that is, we 



r^"^ (p, -P, q, 

set the loop momentum q to zero in the 3 and 
4-point functions. (In this appendix we do not 
indicate explicitly the dependence on (f> of all n- 
point functions in order to alleviate the nota- 
tion.) By doing so, one obtains a closed equa- 
tion for the 2-point function F^ ' {p) , which is the 
object calculated with optimum accuracy at the 
level s — 2. As explained in the main part of 
the article, the general strategy to obtain the 3 
and 4-point functions with comparable accuracy 
is to consider higher orders (s > 2) in the ap- 
proximation scheme. However, in this appendix 
we show that one can already improve the ac- 
curacy of fJ^^' and f[,'^^ simply by exploiting the 

(2) 

information available on F^ (p) . 

Let us consider first the function F^,'^-'. We 
know that, at one-loop and in vanishing fields, 
it has the following structure 



ri.'*^'"°°''(pi,P2,P3,P4) 



f{pi +P2) + f{pi +P3) + f{Pi +P4), 



(Bl) 



where the function f{p) is easily found to be 

fip)^lri'\p,~p,o,o)-lri'\oAO,o). 



(B2) 
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Since T'j^\p, -p, 0,0) = dlTf>(p) (for constant 
field 0), we arrive at the following expression for 
the 4-point function in terms of the 2-point func- 
tion V^^\p): 



k 



1 



P3) 



■dlrfip,. 



P2) 



Pi) 



1 



'2^0rf (0). 



(B3) 



Note that this expression is, by construction, 
symmetric under the exchange of the external 
legs, and it is 1-loop exact at zero external field. 

For the function , one can extract the fol- 
lowing equivalent in the limit of vanishing field: 

d^Tf{p,q,l) 
T^^\p,q,l,0) 



6=0 



= 



(B4) 



(4) 



Then, by using the approximation above for Fj^, 
(IB3p one obtains the following expression for 
F^'^'', whose zero field equivalent is exact at one- 
loop: 



+^9,Ff (0-^a,Ff (0). 



(B5) 



At this point, we note that we may use the new 
expressions that we have obtained for F^. ' and 



F^'*^ in the flow equation for F^^''. Since these 

n-point functions are now one-loop exact, the re- 

(2) 

suiting approximation for FJ^, will be 2-loop ex- 
act in zero external field. This yields therefore 
an improvement of the BMW approximation, in 
particular in the high momentum region where 
we know that it loses accuracy. 

Consider then Eq. (fTTj) for F^ , and re- write it 
in terms of Afc(p): 

5tAfe(p;p)= / dtRk{q)Gl(q) 

Jq 

x[[T['\p,q,-p-q)]^Gk{q+p) 
-[r^^\0, q,-q)]'Gk{q) 



(2) 



-p,q,-q)-Ti^\0,0,q,~q)]}. 



(B6) 



Next, perform the substitutions (jB5l) and 



ri'Hp, -P, q, -q) ^ ^dlTl^'ip + q,~P~ q) 



1 



,2t.(2). 



+ ^dlr'i^\p-q,-p + q). (B7) 



One then gets 



dtAkip) = 2pHip) - -Lip), (B8) 



with 



H{p)^ j dtRk{q)Gl{q)[Gk{q+p)[\K{p) 

+ \K{q) + ^AUp + q)+ + 2pT4"f 
-Gfc(g) [AUg) + Wi; + 2pVi!f], (B9) 



and 

L{p)^ I dtRk{q)Gl{q)U'^{p + q) + 

2pAUp + q)- A',{q) - 2pAt{q)}. (BIO) 

It is not difficult to generalize these expressions to 
the 0{N) model with arbitrary N. However, we 
do not present these here because, in spite of the 
good properties presented above, this extended 
version of the BMW approximation proves to be 
numerically unstable and we have not been able 
to solve the corresponding equations with sim- 
ple techniques. A further analysis, using more 
elaborate numerical techniques, is called for. 



Appendix C: Integrals 

In this appendix, we give details on the calcu- 
lation of the integrals /„(fc; p) and Jn{p, k; p). 

In the case of the integral /„(fc; p), since Gk{q) 
(in a uniform external field) depends only on q"^, 
the angular integral is straightforward. One gets 



where 



(27r)'^ Jo 



dqq''-' dtRkiq) Gl{q-p), 



Sd — 



F(d/2)' 



d{2Tr)d' 



(CI) 



(C2) 



In the case of the integral Jn{p', k; p), the presence 
of the external momentum p makes the angular 
integral more involved. 



1. Angular integrations 

Consider integrals generically of the form 

-l{q)Fi\p + q\)=Iip). (C3) 
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One can proceed as follows 
Ap) = I 9{q)F{\p + q\) 



Sd-i 



i^^r Jo 
r 

X d9 sin' 
Jo 

Sd-i 



d-2 



9 F[\J-p^ + + 2pqcos( 



dq g 



d-1 



9{q) 



{•^T^Y Jo " P 

rP+Q 

X / dUJdi^,P,q)F{0, (C4) 

J\p-q\ 

where we made the change of variables ^ = 
\/p'^ + q"^ + 2pqcos9, and 



Jd{^,P,q) 



e~p' 



2pq 



(C5) 



The interest of the last formula (|C4[) lies in the 
fact that the needed integration points belong to 
the grid so that the integral can be calculated nu- 
merically without the need of interpolation. Fur- 
thermore, this method is particularly convenient 
in d = 3 because the Jacobian (jCSp is then triv- 
ial. 



3. Small momentum limits 

The integral J„(p;fc;p) is regular when p — )• 
0. However, the expression given by the angular 
integration does not make this manifest. To get 
the small p behavior of the generic integral (IC3[) . 
one can expand directly F{\p + q|) — F{q) in the 
first line of eq. ([C4l) : 

F(|p + q|) - F{q) = ^P-iL±^a,F(q) 



2q 



(P - q)^ 

2q^ ^ ^« 



1 



d',F{q) - -d,F{q) 

q 



(C7) 



Then one can use, with the brackets denoting 
angular averages, 

(/('Z))-^/rff^d/('Z) = /('Z) 
((q-p)/('Z)) = ^y' dfld (q-p),/('7)=0 
((q-p)V(g))-^y'rff^d(q-p)V(9) 



9 9 



fiq), 



(C8) 



to obtain 



2. Dimensions less than- 3 

The Jacobian (jC5| . is unity in d = 3, and reg- 
ular for d > 3, but becomes singular for d < 3. 
More precisely, for d < 3 it diverges when ^ ap- 
proaches the boundaries of its integration domain 

= p + q or \p ^ q\). Even if the integral even- 
tually converges, this divergence is the source of 
numerical difficulties. 

We then use for d < 3 a different strategy, 
based on Cartesian variables. We define qi as the 
component of q along p, and proceed as follows 

X / dgi g{q)FUp-^ + 9? + ?! + 2pqi) , 

(C6) 

with q the modulus of the vector q: g = \/ q f + q2 
and \p + q\ = y/p'^ + qf + q^ + 2pqi. 

This expression has no singularities for d > 2 
but it requires multiple interpolations that make 
the numerics more involved than in d > 3. 



(nip + qD-n?)) 

= (d'^F{q) + ^d.Fiq)^ + 0{p% 

(C9) 

and 

Ap) - 1(0) 

x[d^^Fiq) + ^dgF{q))+0{p') 

= \Kd dqg{q)d,{q''-'dgF{q)j 

(CIO) 

Appendix D: Generalization to 0{N) models 

In this appendix the s — 2 BMW approxima- 
tion and the corresponding flow equations are 
presented for 0{N) models. The exact flow of 
the 2-point function in a constant external field 
reads (we omit the renormalization group param- 
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eter k in this appendix for notational simplicity) : 

Jq 

^^ajhiP^ 9' -P ~ 9, 4>)Ghi{q+p, 0) 
+ 9, -q, <f>)G„in{q, 4>) 

^\Gij (<?, 0)rilj7i b, -p, g, -g, (f))Ghn {q, </>) | , 

(Dl) 

where a,b, . . . denote 0{N) indices and cp a A^- 
component uniform field. Within the BMW ap- 
proximation, we make the substitutions; 



(D2) 



In order to manifestly preserve the 0{N) sym- 
metry along the flow, the regulator ASk has to 
be an 0{N) scalar and, accordingly, the cut-off 
function a tensor 

The symmetry of the theory also implies that 
the matrix of 2-point functions can be written 
in terms of two independent tensors. We chose 
to write it in the form 

(2) 

Kb (P^ ~P^ 0) (p, p)5ab + (l)a(j)bTB {P, p), 

(D3) 

with p — J2a 4>a4>a- This form turns out to be 
convenient in the limit p — 0. 

The symmetry also allows us to write the prop- 
agator in this equation in terms of its longitudi- 
nal and transverse components with respect to 
the external field 



Gab{p,(t>)^GT{p',p) Sab 



<t>a4'b 

2p 



+Gl{p',p)- 



2p ■ 



(D4) 



It is easy to find the relationship between these 
propagators and and Tb 

Gt\p,p)^Ta{p,P) + R[p), (D5) 
Gl\p,p) =Ta{p,p)+2pTb{p,p) + R{p). (D6) 

Using the definition (jD3|) of the functions and 
Tb , as well as the form given above for the prop- 
agators, one can decompose the fiow equation 
(IDip in two equations for F^i and F^. 



As in the case TV = 1, we introduce the func- 
tions 

AAip,p) =rAip,p)-p^ -TAip = 0,p), (D7) 
Ab(p,p)=Fb(p,p)-Fb(p = 0,p). (D8) 

Notice that at bare level Ta{p,p) — Ta{p = 
0,p) = p2 while Fb(p,p) - Fb(p = 0,p) = 0, 
which explains the difference between the two 
definitions. In terms of these functions, F^ and 
Tb read 

TAip,p)=P^ + AAip,p) + V\ (D9) 
Tb{p,p)=Ab{p,p) + V'\ (DIO) 

where the primes denote derivatives with respect 
to p. The equations for A^i and read: 

dtAA{p,p)^2p{4^{A'A + V"f 

+ j^iAB + v"f ~ {i^'^ + irw"^} 



li'^{A'^ + 2pA'^) 
iriiN-l)A'^ + 2AB), 



(Dll) 



dtABip, p) = Ji ' (N - 1)(Ab + V"r 

- J^^iA'j, + V"f - J^i^B + V"f 
+ 4^[{A'A + 2AB + W'f 
+ 4p(A's + V") {A'a + 2Ab + 31/") 
+ Ap\A's+V"'f}-\ir{N -l)A's 



1 



J^^{bA's + 2pA'i,) 

^liN-i)ir-h'^-ir)v"' 

-I^'^iSV" + 2pV"'f + AbIa 

(D12) 

where we have omitted the p and p dependences 
on the right hand side for compactness. We have 
introduced the integrals {n > 1) 

J:^{p,p)^ I dtR{q)G"^-'{q,p)Gp{p + q,p), 

lf{p)^jfip^O,p), (D13) 

with a, j3 standing either for L (longitudinal) or 
T (transversal). For n = 1 we set 

/i = (7V-l)/r(p)+/i"(p)- (D14) 
It turns out to be useful to also introduce the 
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integral 



Ia{p) = I dtR{q){GL{q,p)+GT{q,p)) 

X GL{q,p)GT{q,p), 



and, in intermediate steps, we have used the iden- 
tity 



■{GUq, P) - Gl{q, p)) = 2GL{q, p)GT{q, p) 
xrB(g,p)(GL(g,p)+GT(g,p)), (D15) 



which allows us to handle expressions that are 
manifestly regular for p = 0. 

As said in the main text, an accurate study of 
the critical regime requires to use dimensionless 



variables. Using again W{p) = V'{p) we define: 
k'Zuif +KA{prp)) =p^ + Aa{p,p), (D16) 



(D17) 



We also have to use the dimensionless functions 
corresponding to (jD13l) : 



jf{prp)^jf{p.pf^^' 



(D18) 



if{p) = lf{p) 



For numerical reasons, as explained in the 
main text for the = 1 case, we study the flow 
of 

YA{prp)^^, Yb{p,p)^^. (D19) 

pZ pZ 

The dimensionless flow equations can then be 
calculated from equations (jDll[) and (|D12[) : 



dtYAip, p) = v{Ya + 1) + P-g^ + {d-2 + i^Ypy:, 

- + '^PYa) - {{N - + 2Yb) , (D20) 

dtYB iP, p) = (d - 2 + 2r^)YB + p-^ + (d - 2 + Ti)pY's 

+ {N - i)jrf{YB + + Js''{f{n + 2Yb + '-fr 

- jf v(^; + Ji^P^n + -^ybw' + ^) 
- 1^ + 2pw^")' - ((^ - i)^r - - ir) 

- - - ^^"(5^ + 2pY^) + YbIa, (D21) 



with the primes now denoting derivatives w.r.t. The flow equation for the potential, which 

p, and we have omitted the p and p dependences 
on the right hand side for compactness. 
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reads 



We also define 



dtV{p) = -h{p), 



(D22) 



allows us to derive an equation for the dimen- 
sionless derivative of the potential 

dtW{p) = -(2-7y)iy + (d-2 + ,7)pVK' + i^^. 

2 op 

(D23) 

The flow of 77fc follows from fixing a renormal- 
ization condition analogous to Eq. We im- 

pose for all values of fc, 



Ya{po,Po) 0. 



(D24) 



rM = ~vfr{q) - fdgr{q), (D27) 



and the dimensionless propagators 
1 



Gt{p,p) 



Gl{p,p) 



p^{YA + l + r{p)) + W 
1 



(D28) 



+ 1 + 2pYB + r(p)) + W + 2pW' 
(D29) 



from which follow the expressions 

/•OO 

ifip) = d / dqq''-'n{q)Gl~\q)Gp{q), 
Jo 

(D30) 



The simplest choice is po = and po = 0. It 

leads to 



m = :,{NY\{0, 0) + 2YB{0,0))ir{p = 0), 



(D25) 

where we have used ^^{p = 0) = l2^{p = 0). 

In the case of a generic renormalization point, 
the equation for 77^ is more cumbersome 



+ 2po 



1 + PoY^ 



X iPo^ + (d-2)pof; 



dp 





jLT 

^3 ~1 
PO 


~ 2 . 
PO 




-'3 

Po 


~ 2 ' 

Po 



2po J3^^(pV>"J+2rBVF' 

(D26) 

with all functions evaluated at p = po , p — po- 



■p+q 



dq^rMGliq) 



(D31) 

with Jd{Cp,q) as defined in Eq. (jC5|) . We also 
need the functions 

^^(P)^l dqq^-'[rt{q){GL{q) 

+ GT(g))Gi(g)GT(g)(fs(g)g2 + 

(D32) 



rt{q)Ga{q) 



2dKd (27r)'^ JO 

x9jg'^-i9,-G^(g)) , (D33) 



K"^{p) 



that are used in the small momentum region of 
the flow equations (cf. Appendix [Clin the N — 1 
case) . 
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